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CONTAINMENT PROBLEMS FOR PROJECTIONS OF POLYHEDRA 

AND SPECTRAHEDRA 

KAI KELLNER 


Abstract. Spectrahedra are affine sections of the cone of positive semidefinite matrices 
which form a rich class of convex bodies that properly contains that of polyhedra. While 
the class of polyhedra is closed under linear projections, the class of spectrahedra is not. 
In this paper we investigate the problem of deciding containment of projections of poly¬ 
hedra and spectrahedra based on previous works on containment of spectrahedra. The 
main concern is to study these containment problems by formulating them as polyno¬ 
mial nonnegativity problems. This allows to state hierarchies of (sufficient) semidefinite 
conditions by applying (and proving) sophisticated Positivstellensatze. We also extend 
results on a solitary sufficient condition for containment of spectrahedra coming from the 
polyhedral situation as well as connections to the theory of (completely) positive linear 
maps. 


1. Introduction 

A containment problem is the task to decide the set-theoretic inclusion of two given sets. 
In a broader sense this includes, e.g., radii [13] or packing problems [2]. For some classes 
of convex sets there has been strong interest in containment problems. This includes con¬ 
tainment problems of polyhedra and balls [101 and containment of polyhedra na. In re¬ 
cent years, containment problems for spectrahedra, which naturally generalize the class of 
polyhedra, have seen great interest. Ben-Tal and Nemirovski started that investigation by 
developing approximations of uncertain linear matrix inequalities yielding a quantitative 
semidehnite criterion for the so-called matrix cube problem, the decision problem whether 
a cube is contained in a spectrahedron [3]. Helton, Klep, and McCullough studied the ge¬ 
ometry of so-called free spectrahedra (also known as matricial relaxation of spectrahedra) 
including containment problems dZlilHI. They established connections to operator theory, 
namely equivalence between containment of free spectrahedra and positivity of a certain 
linear map. From that they gained a sufficient semidehnite criterion for containment of 
spectrahedra which coincides with the Ben-Tal-Nemirovski criterion when applied to the 
matrix cube problem. Recently Helton, Klep, McCullough, and Schweighofer showed that 
the quantitative version of Ben-Tal-Nemirovski’s criterion is the best possible [T9] . 

The sufficient semidehnite criterion has been reproved by Theobald, Trabandt, and the 
author by a geometric approach. They also provided exact semidehnite characterizations 
for containment in several important cases [23] • In a second work the authors formulated 
the containment problem for spectrahedra as a polynomial nonnegativity question [2B] . 
Based on this formulation they studied a hierarchy of semidehnite programs each serving 
as a sufficient condition for containment coming from Lasserre’s moment approach 1211 and 
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Putinar’s Positivstellensatz |33] . It turned out that the hrst step of the hierarchy is implied 
by the solitary criterion coming from positive linear maps and the geometric approach, 
yielding hnite convergence statements in several cases. In [23] the author considered a 
different but related hierarchy of semidehnite programs based on the Positivstellensatz by 
Hoi and Scherer [23] . As this approach relies on the geometry of the spectrahedra and the 
dehning linear matrix pencils it allows stronger results. Specihcally, all hnite convergence 
results from [26] can be brought forward to this hierarchy and its hrst step coincides 
with the solitary criterion. In addition, using the connection to the theory of positive 
linear maps, hnite convergence is shown for a special family of 2-dimensional bounded 
spectrahedra. 

This paper is concerned with containment problems for projections of polyhedra and 
spectrahedra. Besides the natural question of extending the results for containment of 
polyhedra and spectrahedra to their projections, the paper is motivated by the growing 
attention the geometry of projections of polyhedra and spectrahedra attracted in recent 
years. Among others they have become relevant in many areas like polynomial optimiza¬ 
tion 01121 El], (real) convex algebraic geometry P, |2Q], and extended formulations of 
polytopes [9]. 

Starting point of our considerations are the methods and results discussed above. More 
precisely, we treat possible extensions of the geometric approach, positive linear maps 
and polynomial optimization to the case of projections. The main considerations and 
contributions are the following. 

(1) Although the class of polyhedra is closed under (linear) projections, containment 
problems become more subtle. This is reflected in the fact that the containment 
problem of two projected polyhedra is co-NP-complete (Theorem 2.4). We formu¬ 
late the problem as a bilinear nonnegativity question (Theorem 3.1) and study its 
geometry. 

(2) As the class of spectrahedra is not closed under projections, containment problems 
are even more subtle. However under additional assumptions (which are common 
in semidehnite programming) a similar formulation as in the polyhedral case is 
possible (Theorem 4.1). 

Retreating to the case where only one set is given as a projection, allows to bring forward 
several results from the non-projected case. 

Based on the polyhedral case we deduce a sufficient semidehnite criterion for con- 


(3) 

(4) 


tainment of a projected spectrahedron in a spectrahedron (Theorem 5.5). 

We establish a rehnement of Hol-Scherer’s Positivstellensatz based on the geometry 
of the projected sets (Theorem 5.7). That allows to state a hierarchy of sufficient 
semidehnite conditions to decide containment. The hrst step of the rehned hierar¬ 
chy coincides with the solitary criterion. As a corollary we gain a Positivstellensatz 
for polynomials on projections of polytopes (Proposition |5. 11 ). 

(5) The connection between containment and the concept of positive linear maps can 
be extended to this case (Theorem 5.14). 

The paper is structured as follows. After introducing some relevant notation, we state 
some basics on projections of polyhedra and spectrahedra; see Section [^ We formulate 
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the containment problem for projections of polyhedra as a polynomial nonnegativity ques¬ 
tion in Section and extend it to the case of projections of spectrahedra in Section 
In Section we retreat to the containment problem of a projected spectrahedron in a 
spectrahedron. 


2. Preliminaries 

Let be the space of symmetric k x k matrices with real entries and be the 
closed, convex cone of positive semidehnite k x k matrices. For x = (xi,... ,Xd) denote 
by the space of symmetric k x k matrices with entries in the polynomial ring M[x]. 
For real symmetric matrices AQ,Ai,...,Ad G a linear matrix polynomial A{x) = 
^0 + Ep=i pAp G 5^[x] is called a linear (matrix) pencil. The positivity domain of A[x) 
is dehned as the set of points in for which A{x) is positive semidehnite, 

•S'/! = {x G I A{x) ^ O} , 

where A{x) ^ 0 denotes positive semidehniteness. The closed, convex, and basic closed 
semialgebraic set Sa is called a spectrahedron. 

Every 'H-polyhedron = {x G | a -|- Ax > 0} has a natural representation as a 
spectrahedron called the normal form of the polyhedron Pa as a spectrahedron. 


r k 

'ai(x) 

0 ■ 


X G A(x) = ai(x) = 

[ i=i 

0 

ak(x)_ 

. 0 } 


where aj(x) = (a -f Ax)i for i G [k]. However the converse is not true, i.e., there exist 
nondiagonal pencils describing polyhedra. Deciding whether a given spectrahedron is a 
polyhedron, the so-called Polyhedrality Recognition Problem (PRP), is NP-hard [35] and 
can be reduced to an l-L-in-S containment problem |1|. Note that the normal form used 
here does not coincide with the normal form used in [251 EH] as we do not require the 
constant term a to be the all-ones vector. 

Following the common notation for bounded polyhedra (polytopes), we call a bounded 
spectrahedron a spectratope. 

Denote by vr : —)■ (x, y) ^ x the linear coordinate projection map. By Fourier- 

Motzkin elimination, given an P-polyhedron P = {{x,y) G | a -|- Ax -|- A'y > 0}, 

the projection of P onto the x coordinates is again an "H-polyhedron. Unfortunately, a 
quantiher-free H-description of vr(P) can be exponential in the input size (d, m, k), where 
k is the number of rows in A and H'; see [38] Sections 1.2 and 1.3] and the references 
therein. 

Given a linear pencil A(x, y) G S^[x, y] with x = (xi,..., Xd) and y = (yi,..., ym) for 
some nonnegative integer m, a projection of the spectrahedron Sa is its image under an 
affine map. By an elementary observation, without loss of generality, we can assume that 
the affine projection is a coordinate projection. 

Proposition 2.1 ([TT] Section 2]). If a setT C is the image of a spectrahedron S under 
an affine map, then there exists a linear pencil A(x, y) G S^[x, y] with x = (xi,..., Xd) and 
y = (yi,..., ym) for some nonnegative integer m such that T is a coordinate projection of 
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Sa C Furthermore, if T and S have nonempty interior, then this can be assumed 

for Sa too. 


Due to the proposition, we always assume that the projection of spectrahedron Sa as 
given by the linear pencil A{x, y) G S’^[x, y] with x = (xi,..., Xd) and y = {yi ,..., ?/„,), 
m > 0, is the set 


( 2 , 2 ) 


= { 


X G 


3y G 


A{x,y) F 0} , 


where vr : —>■ denotes the coordinate projection. 

While projections of polyhedra are again polyhedral, this is not true for spectrahedra 
(see, e.g., [6l Section 6.3.1]). Moreover, whereas spectrahedra are basic closed semialge- 
braic sets (the semialgebraic constraints are given by the nonnegativity condition on the 
principal minors), projected spectrahedra are generally not. Though they are semialge¬ 
braic, they are not (basic) closed in general; see Example 4.4 
We state an easy observation for completeness. 


Lemma 2.2. Let A{x,y) G S’^[x,y] be a linear pencil. 

(1) Sa^^ ^ t^{Sa) ^ 0 . 

(2) If Sa is bounded, then ti^Sa) is bounded. 


The converse of part (2) in the previous lemma is not true in general. 

Throughout the paper we use the following notation. The class of projections of "H- 
polyhedra (resp. spectrahedra) is denoted by ttH (resp. tiS). For integers m,n eZ with 
m < n we write [m, n] = {m, m + 1,... ,n} and [m] = [1,..., mj. 


2.1. Complexity of Containment Problems. The computational complexity of con¬ 
tainment problems concerning polyhedra is well-known [TOl [HI [151 CB]- Recently this 
has been extended to spectrahedra [31 [25]. We shortly classify the complexity of several 
containment problems for projections of "H-polytopes and spectrahedra. For more details 
on the complexity classification see [24] . 

Our model of computation is the binary Turing machine: projections of polytopes 
are presented by certain rational numbers, and the size of the input is dehned as the 
length of the binary encoding of the input data (see, e.g., US]). Consider the linear 
projection map tt : —)■ (x, y) h->• x. An Tr^f-polytope 7r(P) is given by a tuple 

[d] m] k] A] A'] a) with d,m,k G N, matrices A G and A' G and a E 

such that vr(P) = {x G | 3?/ G MF : a + Ax + A'y > 0} is bounded. For algorithmic 
questions, a linear pencil is given by a tuple (d; m; fc; Aq, ..., Ad, Aj,..., A(„) with d,m,k E 
N and Aq, ..., Ad, ..., A'^ G rational symmetric matrices such that the projected 
spectrahedron is given by vr(S') = {x G | 3?/ G : A{x,y) P 0}. 

Containment questions for spectrahedra are connected to feasibility questions of semi- 
definite programs in a natural way. A Semidefinite Feasibility Problem (SDFP) is defined 
as the following decision problem (see, e.g., [271 El]). 

Given d,k eN and rational symmetric k x /c-matrices Aq, Ai, ..., Ad, 

decide whether there exists x G such that A(x) P 0. 


(2.3) 
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Equivalently, one can ask whether the spectrahedron Sa is nonempty. Although checking 
positive semidefiniteness can be done in polynomial time by computing a Cholesky fac¬ 
torization, the complexity classification of the problem SDFP is one of the major open 
complexity questions related to semidefinite programming (see [HUM]). Using semidefinite 
programming techniques, a SDFP can be solved efficiently in practice. In our model of 
computation, the binary Turing machine, SDFP is known to be feasible in polynomial 
time if the number of variables d or the matrix size k is fixed m Theorem 7]. 

We first discuss the complexity classification concerning only projections of polytopes. 

Theorem 2.3. Deciding whether a projected 'H-polytope is contained in an 'H-polytope 
can be done in polynomial time. 

Proof. Let vr(P) = {x G | 3i/ G MD : a + Ax + A'y > 0} be a projected P-polytope and 
let Q = {x \ b + Bx > 0} be an P-polytope. Embed Q into by Q' = {(x, y) G 
]^d+m I -|- Oy > 0}. Then the containment problem n{P) C Q is equivalent to the 

'H-'m-'H containment problem P C Q'. The statement then follows from [T3]. □ 

In the latter theorem, the statement does not differ from the non-projected case. The 
next theorem shows a significant change in the complexity classification when the outer 
set is a projected P-polytopes. 

Theorem 2.4. Deciding whether an (projected) PL-polytope is contained in a projected 
PL-polytope is co-NP-complete. 

Proof. Consider a V-polytope. It has a representation as the projection of an P-polytope 
polynomial in the input data. Thus the containment problem Pi-in-TiPL is co-NP-hard 
since PL-in-V is co-NP-complete. It is also in the class co-NP since given a certificate for 
’P not in TiPij i.e., a point p, one can test whether p E PL and p ^ tiPL by evaluating 
the linear constraints of PL (all have to be satisfied) and by solving a linear feasibility 
problem which both is in P by |3Sl Theorem 13.4]. Therefore PL-'m.-TcPL is co-NP-complete. 
Obviously, the proof remains valid when passing to vr'H-in-Tr'H. □ 

In the remaining part, we study the complexity of containment problems involving 
projections of spectrahedra. 

As the complexity of SDFP is unknown, the subsequent statement on containment of a 
spectrahedron in an H-polytope does not give a complete answer concerning polynomial 
solvability of this containment question in the Turing machine model. 

Theorem 2.5. The problem of deciding whether the projection of a spectrahedron is con¬ 
tained in an PL-polytope can be formulated by the complement of semidefinite feasibility 
problems (involving also strict inegualities), whose sizes are polynomial in the description 
size of the input data. 

Proof. Consider a spectrahedron Sa given by the linear matrix pencil A{x,y) and the 
coordinate projection of Sa onto the x-variables ti^Sa)- Given an H-polytope P = {x G 
I 6 -|- Bx >0} with 6 G Q* and B E construct for each i E [/] the SDFP 

{b + Bx)i < 0, A(x, p) P 0 

involving a strict inequality. Then ti{Sa) ^ P if one of the I SDFPs is not solvable. □ 
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While the 7r5-in-'H containment problem is efficiently solvable in practice, the situation 
changes if the outer set is given as the projection of an "H-polytope. 


Theorem 2.6. 

(1) Deciding whether an (projected) 'H-polytope or a (projected) spectrahedron is con¬ 
tained in a (projected) spectrahedron is co-NP-hard. 

(2) Deciding whether a (projected) spectrahedron is contained in the projection of an 
H-polytope is co-NP-hard. 


Proof. Since the problem H-in-S is co-NP-hard (see |3l Proposition 4.1] and [25l Theorem 
3.4]), deciding whether a projected "H-polytope or projected spectrahedron is contained 
in a (projected) spectrahedron is co-NP-hard as well. This is part (1) of the theorem. 

□ 


Parts (2) is a consequences of Theorem 2.4 


2.2. Hol-Scherer’s Positivstellensatz. Consider a symmetric matrix polynomial G = 
G{x) G iS^[a:] in the variables x = (xi,... ,Xd), i.e., a symmetric matrix whose entries he 
in the polynomial ring M[x]. We say G has degree t if the maximum degree of the entries 
is t, i.e., t = max{deg(Gij) \ i,j E [fc]}. 

For matrices M = G and N G 5^, dehne 

i 

(2.4) (M, N), := ({M,„ iV))'^ E,j ■ {M.„ N ), 


where Eij denotes the I x /-matrix with one in the (f,j)th entry and zero otherwise, 
and {•, •) is the Euclidean inner product for matrices. We refer to (2.4) as the Ith scalar 
product. It can be seen as a generalization of the Gram matrix representation of a positive 
semidehnite matrix. Indeed, for positive semidehnite matrices M and N the I x /-matrix 
{M,N)i is positive semidehnite as well [22] . 

For any positive integer /, dehne the quadratic module generated by G{x) 

(2.5) M\G) = {So{x) + (S(x),G(x)}i | So(^) e S'[x], S(x) G S^'[x]} C 5'[x], 

where S^[x] C is the set of sum of squares k x /^-matrix polynomials. A matrix 
polynomial S = S{x) G iS^[x] is called sum of sguares (sos-matrix for short) if it has 
a decomposition S = U{x)U{xY with U{x) G M^^™[x] for some positive integer m. 
Equivalently, S has the form (H ® [x\t)'^Z{H ® [x]*), where [x\t denotes the monomial 
basis in X up to / = max{deg(S'jj(x))/2 | i, j G [/c]} and Z is a positive semidehnite 
matrix of appropriate size. For k = 1, S is called a sos-polynomial. Checking whether a 
matrix polynomial is a sos-matrix is an SDFP (2.3). Obviously, every element in M^(G) 
is positive semidehnite on the semialgebraic set Sg ■= {x G | G(x) ^ 0}. Hoi and 
Scherer [23] showed that for matrix polynomials positive dehnite on Sg the converse is 
true under the Archimedeanness condition. 

We state the desired Positivstellensatz of Hoi and Scherer. See [28| for an alternative 
proof by Klep and Schweighofer using the concept of pure states. 
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Proposition 2.7 ([23l Corollary 1]). Let I be a positive integer and let Sq = {x G 
I G{x) y 0} for a matrix polynomial G G 5^[x]. If the quadratic module M^(G) is 
Archimedean, then it contains every matrix polynomial F G positive definite on Sq- 

By restricting to diagonal matrix polynomials G and / = 1, one gets the Positivstellen- 
satz of Putinar [33] as a corollary. 

Corollary 2.8. Let G = {gi,... ,gk} F M[x] and Sq = {x & \ g > 0 ^g & G}. If the 

quadratic module 

M{G) = < so{x) + '^Si{x)gi{x) \ so,si ,..., Sm e T.^[x] 

[ i=i 

is Archimedean, then it contains every polynomial / G M[x] positive on Sq- 

Interestingly, the usual quadratic module M.^{G) is Archimedean if and only if AA\G) 
is for any positive integer 1. 

Proposition 2.9. The following two statements are equivalent. 

(1) For some positive integer I, the quadratic module A4^(G) is Archimedean. 

(2) For all positive integers I, the quadratic module M.\G) is Archimedean. 
Furthermore, assume G is a linear pencil. Then AA\G) for any positive integer I is 
Archimedean if and only if the spectrahedron So is hounded. 

The equivalence of the hrst two statements was proved by Helton, Klep, and McCul¬ 
lough for monic linear matrix pencils in the language of their matricial relaxation; see [T8l 
Lemma 6.9]. We recapitulate the proof and extend it to quadratic modules generated by 
arbitrary matrix polynomials. 

Proof. The implication (2) (1) is obvious. To show the reverse implication, note hrst 

that Ai\G) is Archimedean if and only if [N — x'^x)Ii G AA\G) for some positive integer 
N. Let m G N be arbitrary but hxed. We have to show that {N — x^x)Im G Ad™'(G). 
Denote by En the m x m-matrix with one in the entry (1,1) and zero elsewhere and let Q 
be the I x m-matrix with one in the entry (1,1) and zero elsewhere. Clearly, Eu = Q^Q. 
Let {N — x'^x)Ii = Sq + {S, G)i with S = be the desired sos-representation. 

Setting §0 := Q^SqQ = G and S = Eu ® An G TF[x\, we get 

{N-x'^x)E^r = Q^{N-x^x)IiQ = Q^(Ao +(A,G),)Q = Ao+.En (An,G) = Ao + (A,G)™. 

Applying the same to Eu for i G [m] and using additivity of the quadratic module A4"^(G) 
yields {N — x'^x)Im G AA"^[G). 

The last statement follows from [301 Corollary 4.4.2] (see also [23]) together with the 
shown equivalence. □ 


3. A Bilinear Formulation of the tt'H-in-tt'H Containment Problem 
F or a G A G A' G and 6 G M', Be B' G let 

7t{Pa) = {t G I 3i/ G : a + Ax + A'y > O} 

^ and 'x{Pb) = {x G | 3i/' G : h + Bx + B'y' >0} 


KAI KELLNER 


be projections of the "H-polyhedra Pa and Pb, respectively. Note that both t^{Pa) and 
71 {Pb) are P-polyhedra themselves (and thus closed sets). A quantiher-free "H-description 
however can be exponential in the input size {d,m,k) respectively {d,n,l); cf. Section]^ 
Our starting point is the formulation of the containment problem as a bilinear feasi¬ 
bility problem. Interestingly, the projection variables y' of the outer polyhedron do not 
appear in the feasibility system (or the optimization version below) only the corresponding 
coefficients B'. 


Theorem 3.1. Let ti{Pa) o^nd 7t{Pb) be as defined in (3.1) and 7r(P^) be nonempty. 

(1) 77{Pa) is contained in 7 i{Pb) if and only if 

z^{b + Bx) > 0 on 7t{Pa) x (ker(P'^) fl . 

(2) Let ker(P'^) fl = span(P')-‘- fl ^ {0}. Then vr(P^) C 7i{Pb) if and only if 

z^{b -f Bx) > 0 on 7i{Pa) x (ker(P'^) fl A*) , 
where = { 2 ; G | = 1 , z > 0 } is the l-simplex. 

The additional assumption on the kernel of seems to be somewhat artihcial, how¬ 
ever, if the projection of Pb to the ^-coordinates is bounded, then the condition holds. 
The two main advantages of part (2) in Theorem 3.1 are the boundedness of the 2 ; variables 
and that the condition z^{b + Bx) > 0 is indeed an inequality. (Note that in part (1), 
containment is equivalent to z^{h+Bx) = 0 on t^{Pa) x (ker(P''^)nM(^) as (x, z) = (x, 0 ) is 
a feasible solution for all x G t^IPa)-) The next lemma serves as a hrst step in a geometric 
interpretation of this precondition. 


Lemma 3.2. Let 7 t{Pb) be as in (3.1). Then ker(P'^) fl M(,_ = span(P') 
and only z/span(P') fl 7 ^ 0. In this case, 7 i{Pb) = 


n ]R+ = {0} if 


In particular, if 7 i{Pb) is bounded, then ker(P'^) fl = span(P')-‘- fl 7 ^ {0}. 


Proof. The equivalence ker(P'^) n]R(,_ = span(P')-‘- n]R(,_ = {0} span(P') 7 ^ 


r^)nM'+ = span(P')^nM^ 
is easy to see. If so, then there exists y' G M" such that B'y' > 0. Thus, for every 
X G M'^, there exists t > 0 sufficiently large such that 6 -|- Px -|- B'{ty') > 0. This implies 


vr(PB) = 


V^. Thus for bounded tt{Pb) we have ker(P'^) fl 7 ^ {0}. 


□ 


3.1 


Before proving Theorem 

7i{Pb) = 1 ^'^” nor the implication '^7i{Pb) is bounded 


we observe that neither the implication “span(P') 7 ^ 

ker(P'^) n ^ {0}” in 


Lemma 3.2| is an equivalence. Example 3.3 also shows that the precondition in part (2) 
of Theorem |3.1| cannot be dropped. 


Example 3.3. 

(1) Consider the polyhedron 



Pi is a pointed polyhedral cone containing the origin in its interior; see Figure UKa). 
We have span(P') fl 7 ^ 0 and thus the intersection of ker(P'^) = ker(l, 1) and the 
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(a) Pi as defined in Example 


3.3 



Figure 1 


nonnegative real nnmbers is zero-dimensional, i.e., ker(i?'^) fl = {0}. Moreover, in 
this case, the restriction to the 1-simplex as in part (2) of Theorem |3.1 is not possible. 
(2) Consider the polyhedron 


A = 


+ 


1 

-1 


X -|- 


-1 

1 


2/>0 


which is nnbonnded and contains the origin in its interior; see Fignre (B). We have 
span(i?') n = 0 and thus ker(P''^) ft ^ {0}. Indeed, for every f > 0, we have 

yiT 


(0, t) G ker(i?'^) On the other hand, ti{Pb) = shows that the reverse of the other 


(and above mentioned) implications in Lemma 3.2 are not equivalences. 


Proof of Theorem \3.1 . 

(1) : 'k{Pa) 2 '^{Pb) if and only if there exists a point x G 7r(Pyi)\7r(PB), i.e., for 

X G 'x{Pa) there exists no y' G M” with b -|- Bx + B'y' > 0. By Farkas’ Lemma [381 
Proposition 1.7] this is equivalent to the existence of a point G with B' = 0 such 

that z^{h + Bx) < 0 holds. Equivalently, there exists {x,z) G tt{Pa) x (ker(i?''^) fl M(,_) 
such that z'^{b + Bx) < 0. 

(2) : If there exists {x,z) G 'x{Pa) x (ker(P'^) 0 IR+) such that z'^{b Bx) < 0, then 

z ^ 0 and thus z''^(b -|- Bx) < 0 for z' = -j^ > 0 with \z'\ = Yl\=i ^ 

Assume z'^{b + Bx) > 0 holds for all (t, z) G tt{Pa) x (ker(P'^) 0 IP^+). By assumption, 
there exists 0 7 ^ ^ G ker(P'^)nM()_. Applying the same scaling as above yields z'^{b+Bx) > 
0 for every z G A*, implying the claim. □ 

Note that ker(P'^) fl is a polytope and is intrinsically linked to the polar of tt{Pb)- 
Namely, it is the set of convex combinations of the columns in that are equal to the 
origin, i.e., 0 = B'^z with 1 = ifz and z > 0. 


Consider the optimization version of Theorem 3.1 


(3.2) 


inf z'^{b + Bx) 

s.t. {x, y, z) E Pa X (ker(P''^) fl Af') . 
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Assuming nonemptyness of ker(i?'^) fl A^, Theorem 
and only if the infimum is nonnegative. 


3.1 


implies that '^{Pa) ^ '^{Pb) if 


Replacing the nonnegativity constraints in (3.2) by sos constraints results in a hierarchy 


of SDFPs to decide the vr'H-in-vr'H containment problem, 

fi{t) := sup fi 

s.t. z^{h + Bx) — /i G + /, 


(3.3) 


where M^ and / denote the quadratic module generated by the inequality constraints and 
the ideal generated by the equality constraints, respectively. 

Under the assumptions in part (2) of Theorem 3.1, applying Putinar’s Positivstellensatz, 
Corollary 2.8, to problem (3.3|) (also allowing equality constraints), the sequence /r(f) 


convergences asymptotically to the optimal value of (3.2) for t ^ oo 


4. A Bilinear Formulation of the nS-m-nS Containment Problem 


Throughout the section, let 

d m 

A{x, y) = Ao + ^ AiXi + ^ A'-yj e S'"[x, y] 

i=l j=l 

d n 

and B{x, y') = Bo + ^ BiXi + ^ R'u/' G [x, y'] 

i=l j=l 

be linear pencils with y = {yi,... ,ym) and y' = (?/(,...,?/(j) for n > 1. Denote the 
projection of the corresponding spectrahedra onto the a:-variables by 

^{Sa) = {x G I 31/ G M”* : A{x,y) P O} 

and '^{Sb) = {x G | 3y' G M” : B{x,y') P O} . 

Recall from Section that the projection of a spectrahedron is not necessarily closed and 
thus, in general, not a spectrahedron itself. 

Dehne B = span{Rj,..., B'^} and recall the equivalence 

{B', Z) = 0 Vi G [n] ^ Z G B^. 


The TTiS-in-TTiS containment problem is slightly more involved than the vr'H-in-vr'H prob¬ 
lem as the projection of a spectrahedron fails to be closed in general. We state an extension 
of Theorem 3.1 to the vriS-in-vriS containment problem. 


Theorem 4.1. Let A{x,y) G S^[x,y] and B{x,y') G S’‘[x,y'] be linear pencils such that 
7i{Sa) ^ 0. 

(1) 7r{SA) P clvr(S's) if and only if (H(x, 0), Z) > 0 on vr(S'A) x fl iS^). 

(2) Assume that the condition A 0 ~ 0 holds for all y'. Then 

the closure in part (1) can he dropped. 


As in the vrH-in-vrH problem, the projection variables y' of the outer spectrahedron do 
not appear in the polynomial formulation, only the corresponding coefficient matrices. 
We use the following Farkas type lemmas to prove Theorem E3 
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Lemma 4.2 ([371 Theorem 2.22]). Let A{x) G 5^[a:] he a linear pencil and denote by 
A{x) = pure-linear part. Then exactly one of the following two systems has 

a solution. 


(4.1) > 0 e 3x e : ||/lo - 4l'o|| < + A{x) G S’f 

(4.2) 3Z : Z ^ 0, {Ai, Z) = 0 Vi G [d], (/Iq, Z) < 0 


The Farkas type lemmas for cones (and thus the theory of semidehnite programming) 
lack in the fact that the linear image of the cone of positive semidehnite matrices is not 
closed in general. Additional conditions which lead to more clean formulations are called 
constraint qualification. 

Lemma 4.3 ([71 Example 5.14]). Let A{x) G iS^[x] he a linear pencil. Assume 




Xj >- 0 




Xi = 0 


i=l 


2 = 1 


holds for any x. Then either ( |4.2[ ) has a solution or Sa is nonempty. 

If Ai,..., Ad are linearly independent, then the above condition can be replaced by 

^ 0 X = 0. 


Proof of Theorem f.l 


(1) : Assume vr(S'A) C cl7r(S's). Let x G 71(5'^)- Then there exists a sequence {xi,y'f)i C 
Sb such that limj^oo Xi = x. For all Z G S-*- fl it holds that 

(B(x,0),Z) = lim {B{xi,y'i),Z) >0. 

2^CX0 

Since x G 7r(S'^) is arbitrary, {B{x, 0), Z) is nonnegative on 71(5'^) x fl iS]^). 

Assume {B{x,0), Z) > 0 on 7 r(S'yi) x {B^ n5^). Let x G 71 ( 6 *^) be hxed but arbitrary 
and set Bq = B{x,0). By Lemma 4.2, there exist Bq G and y' G M" such that 
Bq + YAi=i ^ ‘^+ \\Bq — BqW < e for all 5 > 0. By letting £ tend to zero, there 
exists a sequence {y'Ae ^ such that B{x,y'A ^0. As x G 71 ( 6 '^) is arbitrary, 

the claim follows. 

(2) : Assume (i?(x,0),Z) > 0 on 7i{Sa) x {B^ fl 5]^). Let x G ti^Sa) be hxed but 

arbitrary. By Lemma 4.3, the spectrahedron {y' G M"' | Bq + YAi=i ^iVi — 0} nonempty. 
Thus there exists y' G M"' such that B{x,y') ^0. □ 


Unfortunately, the if-part in Theorem 4.1 (1) without taking the closure is generally 
not true as the next example shows. 


Example 4.4. Consider the linear pencil 


-y'l X 

0 


■Q 

X 

0 ■ 


■-1 

0 

o' 


'o 

0 

o' 

X 1 - 7/2 

0 

= 

X 

1 

0 

+ y'i 

0 

0 

0 

+ 2/2 

0 

-1 

0 

0 0 

-X + 7/2_ 


0 

0 

—X 


0 

0 

0 


0 

0 

1 


and let A(x) be the univariate linear pencil 


A(x) 


1 — X 0 


'1 

o' 


■-1 O' 

0 1 + X 

— 

0 

1 

+ X 

0 1 
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describing the interval Sa = !]• By inspecting the principal minors of B, the spec- 

trahedron Sb has the form {{x,y') G | < 0, x < 1/2 < 1, y[{l — y' 2 ) + < 0}. For 

X = 1, the second condition implies 1/2 = 1 and thus the third condition reads as x^ < 0, 
a contradiction. Thus Sa ^ t^{,Sb) = (— 00 ,1). 

For every Z it holds that 

0 = (Z, B'^) = —Z\\ Z\2 = 0, 0 = (Z, B'^ = Z 33 — Z 22 


implying (i?(x, 0 ), Z) = Z 22 + x(-Z 33 + 2 Z 12 ) = Z 22 (l - x) > 0 for all x G Sa- 

It should not be surprising that the constraint qualihcation on the pencil B (x, y') is not 
satished. Indeed, for {y[,y 2 ) = (//i^O) with y[ < 0, 


B[y[ + B'-^y'^ 


'-y'l 

X 

0 


'-y'l 

0 

o' 

X 

-1/2 

0 

= 

0 

0 

0 

0 

0 

-X + y' 2 _ 


0 

0 

0 


is positive semidehnite but not identically zero. 


is the unboundedness of 
, can be replaced by 

the spectrahedral analog of the simplex. 


An issue when considering the practical utility of Theorem 
the set B^ n iSi. Under an analog condition as in Theorem 


4.1 


3.1 


Corollary 4.5. Let A{x,y) G S'^[x,y] and B{x,y') G S’‘[x,y'] be linear pencils such that 
7 i{Sa) 7 ^ 0- Assume B-^ {0}. Then ttISa) ^ cl 7 r(S'B) if and only if 

(5(x, 0), Z) > 0 on 7r(^A) x {B^ n T^) , 


where = {Z G 5^^ | (J^, Z) = 1} zs the l-spectraplex. 


Proof. Since B^ B^ H the “only if”-part follows from Theorem 4.1 

For the converse, hrst suppose there exists (x, Z) G t^{Sa) x {B^ fl S^) such that 
(i?(x,0),Z) < 0. Then 0 7 ^ Z G and thus tr(Z) = {h, Z) > 0. This implies 

(S(x,0),Z') < 0 for Z' = ^ with tr(Z') = (A,Z') = ^(/^,Z) = 1. 

Assume (i?(x, 0),Z) > 0 on 7i{Sa) x n5^). By assumption, there exists 0 7 ^ Z G 
B^ n n 5^. Applying the above scaling, the claim follows. □ 


We state an analogue to Lemma 3.2 


As the proof is very similar, we skip it here. 


Lemma 4.6. We have = {0} if and only 7 ^ 0. In this case, 7 r(S'B) = M'^. 

In particular, ifn^SB) is bounded, then B^ fl 7 ^ {0}. 


Restricting the 7 r 5 -in- 7 r 5 containment problem to the special case irS-in-TTl-L allows to 


state improved versions of Theorem 4.1 and Corollary 4.5 


Proposition 4.7. Let t^{Pb) be as in (3.1) and let A{x,y) G S’‘[x,y] be a linear pencil. 

(1) tt^Sa) C '^{Pb) if and only if z^{b + Bx) > 0 on 7 r(S'yi) x (ker(i?''^) 

(2) Assume ker(R''^) nM(,_ 7 ^ {0}. Then ^{Sa) ^ ^{Pb) if and only if z^{h + Bx) > 0 
on 7r(S'A) X (ker(i?'’^) n Bf). 
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Proof. tt{Sa) ^ '^{Pb) if and only if there exists x G 7r(S'^) such that fiy' G M"" : h + 
Bp + B'y' > 0. By Farkas’ Lemma [38l Proposition 1.7] this is equivalent to the existence 
with z'^B' = 0 and + Bp) < 0. The claims follow as in the proofs of 

^ □ 


of a z G 
Theorem 4T and Corollary |4.5[ 

We close with an example. 


Example 4.8. Let M be the convex hull of the shifted unit disks dehned by the identities 
1 — (xi + 1)^ — X 2 = 0 and 1 — (xi — 1)^ — x^ = 0, respectively. M is the projection of a 
spectrahedron. Indeed, considering only the hrst disk and shifting it along the segment 
[—1,1] X {0} yields M = {x G | 3i/ G M : 1 — {x — y)'^ — x^ > 0, —l<y<l}. It is 
the projection of the 3-dimensionaI cylinder, see Figure dehned by the linear pencil 


A{x,y) 


I-X2 xi-y 
xi-y I + X 2 


1-y 0 

0 1 + y 



Figure 2. The determinantal variety of A{x,y) with Sa being the grey 


cylinder in the middle of the picture; see Example 4.^ 


The so-called TV screen (see, e.g., [6l Section 6.3.1]) is the projection of the spectrahe¬ 
dron 


Sb = |(j:,2/) e I A{x,y) 
onto the x variables; see Figure 


1 + 1/1 

y2 

© 

1 Xi 

© 

'l 

X2 

+ o| 

y2 

1 -?/i. 


jTi yi 


_X2 

y2_ 

" i 



Figure 3. 


The TV screen as stated in Example 4.8 
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Both 7 r(S'^) and 7i{Sb) are closed but not spectrahedra. We have 

B"*" n T® = I Z G 5^ I ^ ^ Zii = 1, Z 22 = Z\\ + Z 44 , Zge = —‘IZ 


i=l 


For all Z G B-*“ fl T® the objective (i?(a;,0), Z) in Corollary 4.5 can then be written as 

{B{x, 0), Z) = 1 — Z 44 — Zqq + 2 X 1^34 + 2X2^56- 

We want to hnd a pair (x*, Z*) G 7i{Sa) x (B-*- fl T®) such that {B{x, 0),Z) < 0. To this 
end, consider x* = (1 + e, 0) G vr(5yi) for all e G [0,1] and 


Z* = 


0 

o' 


■ 1 

1" 

0 

1 

© 



3. 


. 3 

3 . 


O.SX.S G iS, 


Then {B{x*, 0), Z*) = —< 0 for all e G (0,1]. Thus 7 r(S'^) ^ vr(S's). 

Now interchange the roles of Sa and Sb, i-e., h^Sa) is the TV screen and 7i{Sb) is the 
convex hull of two disks. Then B-*“ fl is the set 


nr = \zeSl\ Y,Zr. 


— 1, 2 Zi2 — Z44 — Z: 


44 — 4^33 


2 = 1 


For all Z G B-*- fl (i?(x, 0), Z) has the form 

(B(x, 0), Z) = (1 — X 2 )Zii + (1 + X 2 )Z 22 + (1 — Xi)Z 33 + (1 + Xi)Z44. 

As l±Xj > 0 for all x G vr(S' 2 i), we have (i?(x, 0), Z) >0 for all (x, Z) G 71 ( 6 '^) x (B-*- flT'^). 
Thus 7r(S'^) C 7t[Sb)- 


5. Sum of Squares Certificates for the ttiS-in-^S Containment Problem 


Retreating to the cases n'H-'m-'H and vTiS-in-iS allows to bring forward several results 
from the non-projected case. We start with the polyhedral situation in Theorem |5.1[ It 


also serves as an algorithmic proof of Theorem 2.3 Afterwards, we state and prove a 


sophisticated Positivstellensatz for the second problem. 


5.1. Prom the 'K'H-in-'H to the TTiS-in-iS Containment Problem. As the proofs of 
the statements in this section are similar to the ones given in [ 2 S], we only stress the 
emerging differences in the proofs. 

Even in the non-projected case, i.e., m = 0, Theorem 5T below is a slight extension of 
a statement in [25]. Namely, here we drop the conditions a = 1^, and 6 = as well as 
the boundedness condition. 


Theorem 5.1. Consider the polyhedra Pa = {{x,y) G | a + Ax + A'y > 0} 7^ 0 and 

Fb = {x G I 6 -I- Bx > 0}. 

(1) 'x{Pa) ^ Pb if o,nd only if there exists a nonnegative matrix C G and a 

nonnegative vector cq G with b = cq + Ca, B = CA, and 0 = CA'. 

(2) Let Pa he a polytope that is not a singleton. Then '^{Pa) ^ Pb if o,nd only if there 
exists a nonnegative matrix C G ^ with b = Ca, B = CA, and 0 = CA!. 
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Testing whether Pa is a singleton is easy as one has to check that the system of equalities 
a + Ax = 0 has a single solution. Certainly, in this situation, checking containment is 
trivial as 7^{Pa) ^ Pb is equivalent to test whether a single point has nonnegative entries. 
The precondition in part (2) of Theorem 5.1, however, cannot be removed in general; see 


part (1) of Example 5.2 


For unbounded polyhedra the additional term cq is required in order for the criterion to 
be exact. Without it, already in the simple case of two half spaces dehned by two parallel 


hyperplanes, the restriction of the condition in part (1) of Theorem 5.1 to part (2) can 


fail to be feasible; see part (2) of Example 5.2 


Example 5.2. 

(1) Consider the polytopes Pa and Pb given by the systems of linear inequalities 



-1 -1 

1 0 

0 1 


a: > 0 and j 2 


+ 


1 0 

-1 -1 

-1 1 


a: > 0, 


respectively. Pa is the singleton {(1,0)} and Pb is a simplex containing Pa- There is 
no matrix C satisfying the conditions in part (2) of Theorem 5.1 (with m = 0). Indeed, 
b = Ca implies 0 = Cu — C 12 and B = CA implies 1 = Bn = (—Cn + C^-, —Cn)., a 


contradiction. A solution to the linear feasibility system in Theorem 5.1 (1) is given by 


Co — Is) C — 


0 1 0 
1 0 0 
1 0 2 


Moreover, it is easy to see that for any Pb = (a; G | 6 + Bx > 0} containing Pa 
containment is certihed if and only if B has the form B = [—b, —b + c] for some vector c. 

(2) Consider the half space given by the linear polynomial a(a;) = 1 — Xi — a; 2 . Let 
b{x) = b + [Bi, B 2 \x be any half space. The condition in part (2) of the Theorem |5.1 


IS 


satished if and only if 6 = c, Bi = —c, B 2 = —c for c > 0. Thus either b{x) = 0 or b{x) 
is a positive multiple of a (a;). 


To prove Theorem 5.1 we use the following affine form of Farkas’ Lemma. 


Lemma 5.3 ( [36l Corollary 7.1h]). Let P = {a; G | a + Ax > 0} be nonempty. 
Then every affine polynomial f G M[a;] nonnegative on P can be written as f{x) = cq + 
Ym=i with nonnegative coefficients Ci- 


Proof of Theorem \5.1\ If P = CA, 0 = CA', and b = Ca (resp. b = cq + Ca) with a 
nonnegative matrix C, for any x G t(P^) we have 

b + Bx + t)y = C {a + Ax + A'y) > 0, 

i.e., ^{Pa) C Pb- 

Conversely, if 'x{Pa) P Pb, then any of the linear polynomials {b + Bx + 0y)i, z G [/], 
is nonnegative on Pa- Hence, by Lemma 5.3, {b + Bx + 0y)i can be written as a linear 
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combination 


k 

{b + Bx + 0y)i = + +Ax + A'y)j 

i=i 


with nonnegative coefficients ch. Comparing coefficients yields bi = c'q + 
i E [/], implying part (1) of the statement. 

To prove the second part, hrst translate both and Pb to the origin. By assumption, 
there exists (x, y) G Pa- Dehne a := a + Ax + A'y and b := b + Bx. Then a > 0 and 
0 e {x G I a + Ax + A'y > 0}, implying 


b = Ca, B = CA, 0 = CA' ^ b = Ca, B = CA, 0 = CA'. 


Thus w.l.o.g. let a > 0. 

Stiemke’s Transposition Theorem [351 Section 7.8] implies the existence of a A > 0 such 
that [A^, A'^jA = 0, and thus 

{a + Ax + A'y) = X^ a = 1 

after an appropriate rescaling. Note that a 7^ 0 as otherwise Pa = {0} is a singleton. By 
multiplying that equation with c'q from above, we obtain nonnegative c'X with Cp (a+ 
Ax + A'y)j = c'q, yielding 

k 

{b + Bx)i = + c''j){a +Ax + A'y)j. 

j=i 


Hence, C = with Cij := ch + c7 is a nonnegative matrix with B = CA, 0 = CA , 

and {Ca)i = = h - c'q + c'q A^a = bi for every i e [1]. □ 

The sufficiency part of Theorem IP can be extended to the case of projected spectra- 
hedra via the normal form (2.1) of a (projected) polyhedron Pa as a (projected) spectra- 
hedron, 

ti^Pa) = {x G I 3i/G M'" : H(x, 1 /) = diag(ai(x, ?/),..., afc(x, 1 /)) ^ O} , 
where ai(x, y) is the ith entry of the vector a + Ax + Ay. 

Corollary 5.4. Let A{x,y) G S^[x,y] and B{x) G 5^[x] be normal forms of polyhe- 
dra (2.1). 

(1) ti^Sa) C Sb if and only if there exist positive semidefinite diagonal matrices C'q, C 
such that 

k k k 

(5.1) i?o = C'q + '^^{AQ)iiCii, Bp = '^^{Ap)iiCii Vp G [d], 0 = '^^^{A^aCu Vp G [mj. 


Z=1 


2=1 


2=1 


(2) Let Sa be a polytope that is not a singleton. '^'{Sa) ^ Sb if and only if system (5.1) 
has a solution with Cq = 0. 
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If the diagonality condition on the matrix C in Corollary |5.41 is dropped, then the above 
SDFP yields a sufficient condition for the vTiS-in-iS containment problem. Subsequently, 
the indeterminate matrix C = is a symmetric kl x A;/-matrix, where the Cij are 

I X /-blocks. 


Theorem 5.5. Let A{x,y) G S^[x,y] and B{x) G be linear peneils. Denote by 

ti{Sa) the coordinate projection of the spectrahedron Sa- If there exist positive semidefinite 
matrices C = G and Cq G such that 

k k k 

(5.2) Bo = Co + '^^^{Ao)ijCij, Bp = {Ap)ijCij Vp G [d], 0 = '^^{A'p)ijCij Vp G [m], 

i,j=l *j=i *j=i 

then 77(5'^) C Sb- 


In the non-projected case, the sufficient semidehnite criterion (5.2) has hrst been devel¬ 


oped by Helton et al. [H] using the theory of positive linear maps (cf. Section 5.4) and 
has been reproofed in [25] by elementary methods. In [23] the author showed that the 
condition is exactly the 0th step of the hierarchy based on truncation of the Hol-Scherer 
quadratic module (2.5). Here we are bringing this forward to the projected case. 


For completeness we state a short proof of Theorem 5.5 based on 


Proof. We have 

d k 

B{x) = Bo+ Y, ^pBp = = C'o + {{A{x, I 

p=l i,j=l 


with I = [//,.. .JiY 


iklxl 


Let X G 7i{Sa)- By dehnition, there exists y G such 


that A{x,y) P 0. Thus the Kronecker product A{x,y) ® C is positive semidehnite. Since 
{{A{x,y))ijCij)’^ is a principal submatrix of A{x,y) (g)C, we have B{x) ^ 0 as well. □ 


Even for the non-projected case, the sufficient semidehnite criterion (5.2) is not neces¬ 
sary for containment in general; see m Section 6.1]. 

5.2. A Sophisticated Positivstellensatz. Consider the linear pencils A{x,y) G 5^[x] 
and B{x) G iS^[x]. Then 7i{Sa) is contained in if and only if B{x) P 0 on 7i{Sa)- If Sa 
is a spectratope, then this is equivalent to B{x) -|- eli G M.\A) for all e > 0, where 

A4'(A) = {^o + (5,A(x,p)), |^oeS'[x,p], SeJ:^%y]} 


is the quadratic module associated to A(x, y) as dehned in (2.5). Clearly, if B{x) G Ai\A) 
for a linear pencil B{x) G 5^[x], then cl7r(S'A) C Sb- Thus truncation of the Ai\A) yields 
a hierarchy of SDFPs to decide 7riS-in-5 containment. 

The drawback of this approach to the vTiS-in-iS containment problem is that it relies on 
the geometry of the spectrahedron Sa rather than its projection, namely the boundedness 
assumption on Sa and the appearance of the projection variables y in the quadratic 
module. In the following, we address this by developing a refinement of Hol-Scherer’s 
Positivstellensatz. Particularly, we can eliminate the variables y in the sense that they 
neither appear in the quadratic module nor in the relaxation. 
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Gouveia and Netzer m derived a Positivstellensatz for polynomials positive on the 
closnre of a projected spectrahedron. 

Proposition 5.6 ([TTl Theorem 5.1]). Let A{x,y) G S^[x,y] be a strictly feasible linear 
pencil. Define the quadratic module 

A^(7rA) = {so + (>5, A(a;, 0)) | {S, A'f) = 0 \/i e [m], sq £ S'gS^[x]}. 

If 'k{Sa) is bounded; then M.{ttA) is Archimedean and contains all polynomials positive 
on the closure o/7r(S'/i). 

Snbsequently, we state and proof an extension to linear pencils positive dehnite on a 
projected spectrahedron. Thereto dehne the qnadratic modnle 

(5.3) Af'(7rA) = {5o +(^,A(a;,0)); | {S, A'fi^ = 0'ii e [m], S-q G S'[x], ^gS^'[x]}. 

It is easy to see that M.\nA) is in fact a qnadratic modnle. Note that AdfinA) does not 
have to be hnitely generated; see m Section 5]. Clearly, every element of AA\ttA) is 
positive semidehnite on the closnre of '^{Sa). 

Theorem 5.7. Let A{x,y) G S^[x,y] be a strictly feasible linear pencil such that ti^Sa) 
is bounded. For I G N the quadratic module AAijrA) is Archimedean and contains every 
matrix polynomial positive definite on cl7r(S'^). 

Proof. By bonndedness of 7r(S'A) there exists iV G N snfficiently large snch that iV ± x* is 
nonnegative on tt^Sa) for all i G [d]. We show that nnder the preconditions in the theorem 
Ai^TiA) contains every linear polynomial nonnegative on 7i{Sa)- Then the qnadratic 
modnle is Archimedean. 

Let b{x) = bo + b^x G M[x]i be a hxed bnt arbitrary affine linear polynomial nonnegative 
on 7i{Sa)- Consider the following primal-dnal pair of SDPs. 

snp {-Ao,Z) 

p* := inf b{x) s.t. {At, Z) = hiMi & [d] 

s.t. A{x; y) Ft) (A', Z) = 0 Vi G [m] 

ZeSl 

Since A(a:, y) is strictly feasible by assnmption, the dnal problem (on the right-hand side) 
has optimal valne p* — bo and attains it; see [H Theorem 2.2]. Since b{x) > 0 on 7i{Sa), 
we have —bo < p* — bo and thns 

bo- Zo = (Ao, Z ), {Ai; Z) =bi\/i e [d], (A', Z) = 0 Vi G [m] 

for some Z E and Zo > 0. Dehne S{x) as the block diagonal kl x /cZ-matrix with I 
copies of Z on its diagonal, i.e., S{x) = and So{x) = zoh- Then 

i 

So{x) + {S{x), A{x, 0)); = zqIi + ^ {Z, A{x, 0))^ = b{x)Ii 

i=i 

and {S{x),A'fii = A') = 0 for i G [m]. This implies b{x) G M^tiA). By Hol- 

Scherer’s Theorem, every matrix polynomial positive dehnite on cl 7i{Sa) is contained in 

M^inA). □ 
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Theorem 5T leads to a refined hierarchy for the 7riS-in-5 containment problem using 
the truncated quadratic module 

(5.4) MliTrA) = {^0 + (^,4l(a:,0)), | = 0 V* G [m], e ^ G Sf [x]} . 


It is evident from the definition of the quadratic modules A4.\A) and AA.\tiA) that the 
latter approach is preferable to the naive way from the theoretical viewpoint (provided 
that A{x,y) is strictly feasible). 


Corollary 5.8. Let A{x,y) G S^[x,y] be a strictly feasible linear pencil such that ti^Sa) 
is bounded and let B{x) G be a linear pencil. 

(1) ti{Sa) C Sb if and only if B(x) + eli G M-^tiA) for alls > 0. 

(2) If B{x) y 0 on 7t{Sa), then B{x) G Ai’‘[7rA). 


The 0-th step of the hierarchy based on (5.4) is exactly the sufficient containment 


criterion stated in Theorem 5.5 Indeed, for f = 0, the constant sos-matrix S equals the 


positive semidefinite matrix C after permuting rows and columns simultaneously. 


Proposition 5.9. Let A{x,y) G S^[x,y] and B{x) G 5*[a:] be linear pencils. Assume 
A{x,y) is strictly feasible. The following are equivalent. 

(1) B{x)eM,{7;A). 

(2) There exist C G and Cq G S’:^ such that 

Bo = c' + (Ao, C), , B, = {A,, C), Vp G [d], 0 = (a;, C')^ Vg G [m]. 

(3) There exist C G S^., Cq G such that 

k k k 

Bo = Cq+ ^^(Ao)pCjj, Bp = ''^^{Ap)ijCij Vp G [d], 0 = '^^{A'^)ijCij Vg G [m]. 


Proof. The equivalence of (1) and (2) follows from the definition of the truncated quadratic 
module by rewriting it as an SDFP. Applying a simultaneous permutation of the rows and 
columns of C in (2) (resp. of C in (3)), the linear systems can easily be transformed. 

For details in the non-projected case see [2ll Theorem 5.1.11] □ 


The proof of Theorem 5.7 evidently yields necessity for the vriS-in-'H containment prob¬ 
lem. This, in particular, shows the (theoretical) effectiveness of the approach based on 
Theorem 15.71 


Theorem 5.10. Let A{x,y) G S^[x,y] be a strictly feasible linear pencil and let the 
coefficients of the linear pencil B{x) G 5^ [a:] be simultaneously congruent to a diagonal 
matrix. 


(1) 'k{Sa) C Sb if and only if B{x) G ^^{ttA). 

(2) Assume Sb is a polytope with nonempty interior. Then 7 i{Sa) ^ Sb if and only if 
B{x) G A4q{7tA) with Bo = 0. 


In particular, the statements (1) and (2) hold for a diagonal linear pencil B{x), i.e., a 
polyhedron in normal form (2.1). 
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In order to prove Theorem 5.10 we use natural adaptions of auxiliary results on the 
behavior of the sufficient containment criterion with regard to block diagonalization and 
transitivity as shown in [25j to the non-projected setting. Using Proposition |5.9| it is easy 
to verify the validity of these statements. 

Proof. As for f = 0 the resulting SDFP is invariant under non-singular congruence trans¬ 
formations of B{x) (see pTl Lemma 5.1.14]), we can retreat to the normal form (2.1) 
B{x) = b‘^{x) E 5^[x] with W{x) = bl + x'^b^ for q E [1]. Denote hj blybl,, b\ the 

coefficients of the linear form b^[x) = {bo + Bx)q. Set 6^ := {bl,... ,bf). 

The proof of Theorem 5T yields certificates 

bl - = (Ao, , (A„ Z^) = bj Vz e [d], (A', Z'') = 0 Vz e [m] 

for some E S’f and Zq > 0. Setting S{x) = 0^=1 Z'^ and S'o(x) = 0g=i Zq, this implies 
part (1) of the statement. 

To prove the second part, let S{x) as before and set S'o(x) to be zero. Then 

i i / 

{S{x),A{x, 0)); = 0 (A(x, 0), Z^)^ = 0 j /o - 4 - 

q=l q=l V 

certifies the containment ^{Sa) ^ S'^/, where B'{x) is defined as 


'^blxi 


i=l 


B'{x) = 0 

q=l 


p=i 


Assuming that S'^ is a polytope, we have Sb' ^ Sb and thus, by transitivity and exactness 
of the initial hierarchy step for polytopes, see Corollary 5A, there is a certificate for the 
containment question tt{Sa) ^ Sb of degree zero with S'o(a;) = 0. □ 


As a special case of Theorem |5.7[ we gain a Positivstellensatz for polynomials on pro¬ 
jected polyhedra having boundedness as its only precondition. 

Proposition 5.11. Let Pa = {{x,y) E | a+Ax+A'y > 0} be a nonempty polyhedron 
such that 'x{Pa) is bounded. The quadratic module 

f k k 

A) = <^ So + X] Si(x)(a -F Ax)i \ ^ Si{x){AA) = 0 Vj G [m], Sq, ..., s^ G S[a;] 


i=l 


i=l 


is Archimedean and contains every polynomial positive on 7 i{Pa). 


Proof. The proof follows from the proof of Theorem |5.7| by retreating to diagonal pencils 
and the fact that strong duality holds for linear programming [36l Corollary 7.1g]. □ 


5.3. Examples. We discuss some academic examples for the hierarchy stated in (5.4). 
All computations are made on a desktop computer with Intel Core 13-2100 @ 3.10 GHz 
and 4 GB of RAM. In the tables, “time” states the time in seconds for setting up the 
problem in YALMIP [32] and solving it with Mosek [1] . 
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Table 1. 


t(S'a) 

tSb 

r 

time 

^(0) 

two disks 

2 -ball 

1.99 

0.7978 

-0.0050 



2 

0.8215 

5.9978 ■ 10-°^ 



2.01 

0.9173 

0.0050 


tSb 




3-ball 

2.23 

0.8690 

-0.0027 



2.2361 

0.7470 

1.4339 ■ 10-°® 



2.24 

0.8803 

0.0018 


Computational test of containment as described in Example 


5.12 


In the examples we consider the optimization version of (5.4) 


fi(t) = sup fl 

s.t. B{x) — fill G 


Letting t tend to inhnity, the sequence of optimal values /r(t) converges to the value 
fi* = sup{/i I B{x) — fill ^ 0 Vx G 71(5'^)} which is nonnegative if and only if ti^Sa) ^ Sb- 
Thus a nonnegative value /r(0) states the existence of a containment certihcate with t = 0. 
A d-dimensional ball is a spectratope Sa given by the linear pencil 


d 

(5.5) A{x) = -^{Ep^d+i + -£'d+i,p) £ 

p=i 

with r > 0. 



Figure 4. 


The convex hull of two disks in a 2-ball as stated in Example 


5.12 


Example 5.12. Consider the convex hull of two disks 71(5^) as dehned in Example 4.8 


and the 2-ball of radius r > 0. It follows from the construction of 7r(S'yi) that it is centrally 
symmetric and that its circumradius is 2; see Figure Up to numerical accuracy, this 
value is computed by our approach; see Table [1} 
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^{Sa) rSs r 

time /r(0) 

TV screen 2-ball 1.18 

0.8514 -0.0078 


2.0564 -1.3621 ■ 10"°® 

1.19 

0.9964 6.6628 ■ 10"°^ 

1.2 

0.9854 0.0090 

tSb 


4-ball 1.55 

1.1036 -0.0024 

W 2 + I 

0.9373 3.1037-lO-o^ 

1.56 

1.0723 0.0040 


Table 2. 


Computational test of containment as described in Example 


5.13 


Example 5.13. Consider the TV screen 7r(S'A) as defined in Example 4.8 Note that 


while the TV screen is centrally symmetric (as its boundary equals the variety dehned by 


the polynomial 1 — x\ — X 2 ), its dehning spectrahedron is not (as the point (1,0,1,0) is 
contained in Sa but its negative (—1, 0, —1, 0) is not). 

As one can see in Table the circumradius of the TV screen is at most ~ 1.1892, 
while the “centrally symmetric circumradius” of Sa is at most a/2 + 1 1.5538. 

Actually, the computed values for the (centrally symmetric) circumradius are exact. 

For p = e Sa we have ||7r(p)||2 = ^ and ||p ||2 = VW+I, implying 

that the circumradius of the TV screen is at least and that \/ a/ 2 + 1 is the smallest 
possible radius of a ball (centered at the origin) containing Sa- 


5.4. Containment of Projected Spectrahedra and Positive Linear Maps. We 

discuss an extension of the connection between positive linear maps and containment of 
spectrahedra (as introduced by Helton et ah [18]; see also [26]) to projected spectrahedra. 

Given two linear pencils A{x,y) G S^[x,y\ and B{x) G iS*[x] with y = (|/i,..., i/^) 
dehne the linear subspaces 

A = span{Ao ,..., A^, A[,..., A'^} and B = span{i?o,..., Bd}. 

Every element in A can be associated to a homogeneous linear pencil A{xo,x,y) G 
S^[xo,x,y] (Ao being the coefficient of xq). The linear pencil 

d m 

A{xq,x, y) := Xo(l © Aq) + Xp{0 © Ap) + yq(0 © A^) 

p=l q=l 


is called the extended linear peneil associated to A(xo, x, y). The associated linear subspace 
is A = span|l © Aq, 0 © Ai,..., 0 © Ad-, 0 © A(,..., 0 © Ajj^j-. 
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For linearly independent Ai,, Ad, A\,..., A'^, let <Fab : ^ be the linear map 
defined by 

© -4o) = Bo, Sab(0 © Ap) = Bp^p e [d], $ab( 0 © Ap = 0 Vp G [m]. 

Since every linear combination 0 = Ao(l © Aq) + X]p=i -^p(0 © ^p) + Ad+g(0 © A'^ 
for real scalars Aq, ..., A^+m yields Aq = 0 , it snffices to assnme the linear independence 
of the coefficient matrices Ai,..., Ad, A'.^,..., A!^ to ensnre that <Fab is well-defined. If 
^ 0 ; Ai,..., Ad, A']^,..., are linearly independent, then we can retreat to the simpler 
map ^ab '■ a ^ B defined by 

^ab{Ap) = Bp yp e [d] and ^AB{^q) = 0 Vg G [m]. 

The next theorem extends the key connection between operator theory and containment 
of spectrahedra to the setting of projections of spectrahedra. 

Theorem 5.14. Let A[x,y) G S^[x,y\ and B{x) G be linear peneils. 

(1) If ^AB or ^ab is positive, then 71 ( 5 '^) © Sb- 

(2) If ti{Sa) 7 ^ 0, then ti{Sa) © Sb implies positivity o/<F. 

(3) If ti{Sa) 7 ^ 0 and Sa is bounded, then 7 i{Sa) © Sb implies positivity o/<h. 

Proof. 

(1) : Let ^ab be positive. For every x G ti^Sa) there exists y G MP snch that A[x, y) © 0, 

i.e., A{x,y) G fl .4,. Then B{x) = B{x) + = ^{A{x,y)) G fl S and 

hence x E Sb- There is no difference in the proof if ^ab is positive. 

( 2 ) : Since the spectrahedra defined by A{x,y) and A{x,y) coincide, their projections 
equal and hence tt {S£) © Sb- Let A{xo,x,y) G (lA. Then xq > 0. 

Case Xo > 0. By scaling the linear pencil with 1/xo the positive semidefiniteness is 
preserved. Thus, l/xoA{xo,x,y) = A{l,x/xo,y/xo) G A A and x/xq G ^{Sa) © Sb- 
Scaling B{x/xo) by Xq yields <h(4(a;o, x, y)) = XqBq + Y.t=i ^pBp = XoB{x/xo) eS^AB. 

Case Xo = 0. If {x,xo) = (0,0), the statement is obvious. Let a; 7 ^ 0. Fix a point 
X E ^{Sa) 7 ^ 0- Then, for some y,y E M™, 4(1, x+tx, y+ty) = 4(1, x, y) +4(0, tx, ty) © 0 
for all t > 0, implying x+tx E ti{Sa) © Sb for all t > 0. Thus x lies in the recession cone of 
^{Sa) which clearly is contained in the recession cone of Sb- Indeed, \B{1, x) + 5(0, x) = 
IB{x + tx) © 0 for all t > 0. By closedness of the cone of positive semidefinite matrices, 
we get 5(0 ,t) © 0. Hence, ^{A{xo,x,y)) = <F(4(0,x,p)) = 5(0 ,t) © 0. 

(3) : Let 4(a;o, x, y) = XqAq + ^p\ + Vq^'q e 5^ n 4. 

Case Xo > 0. This case follows by a similar scaling argument as in part ( 2 ). 

Case Xo < 0. Since ti^Sa) 7 ^ 0, there exists x G 71(5'^) and hence, for some y E M™', 

4(0,a; + \xo\x,y+ \xo\y) © ItqI ■ A{l,x,y) © 0. 

For 4(0, T + \xo\x,y + \xo\y) 7 ^ 0, one has an improving ray of the spectrahedron Sa, in 
contradiction to boundedness of Sa- For A{f), x + \xo\x, y + \xo\y) = 0 , linear independence 
of 4o,.. -,Ad+m implies {x+\xo\x,y+\xo\y) = (0,0). But then a;o4(l, t, p) = A{xo,x,y) © 
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0 together with xq < 0 and A{l,x,y) ^ 0 imply either A{l,x,y) = 0, in contradiction to 
linear independence, or (xo,x) = (0,0). Clearly, in this case, <hyiB(0) = 0. □ 


The linear map ^ab can be represented by an symmetric {k + 1)1 x {k + 1)1 matrix C. 
By expecting the linear eqnations dehning C, it is easy to see that C = Cq + C, where the 
matrix pair Cq, C is from Theorem 5.5 Thns Proposition 5.9 implies the next corollary 
to Theorem 15.141 


Corollary 5.15. Let A{x,y) G S^[x,y] and B{x) G 5^[x] be linear pencils with A{x,y) 
strictly feasible. Then the following are equivalent. 


(1) The map ^ab is completely positive, i.e., 

(2) The solitary criterion (5.2) is feasible. 

(3) B{x) G Ml^inA). 


d + 0. 
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